1 



Status of Lattice Studies of the QCD Phase Diagram 

Owe Philipsen*) 

Institut fiir Theoretische Physik, Westfdlische Wilhelms- Universitdt Miinster, 

- - - 4^149 Miinster, Germany 

00 

' Determining the QCD phase diagram is a pressing task in view of its relevance for nuclear 

' and astro-particle physics programmes. We review the current status of lattice calculations 

, of the phase diagram in the (T, ^_B)-plane for baryon chemical potentials /is < 500 MeV. At 

/is = 0, simulations of staggered fermion actions predict the quark hadron transition to be 
^ ' a crossover in the continuum limit. As a baryon chemical potential is turned on, there is 

I mounting evidence on coarse lattices for the crossover to weaken, rather than turning into a 

true phase transition at a critical point, as predicted by earlier simulations. 
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§1. Introduction 



The QCD phase diagram has been the subject of intense research over the last 
(D \ ten years. Once fully determined, it will locate the regions of different forms of 

7^ , ■ nuclear matter in the parameter space spanned by temperature T and baryon chem- 

ical potential ^b- Based on the fundamental property of asymptotic freedom, one 
^ expects at least three different regions: hadronic (low fiB,T), quark gluon plasma 

■ (high T) and colour-superconducting (high fiB, low T). For chemical potentials ex- 

, ceeding /is > 1 GeV, the situation may be more complicated with possible additional 

^ ' phases.^-* 

Unfortunately, a quantitative calculation of the phase diagram from first prin- 
ciples is extraordinarily difficult. Since QCD is strongly coupled on scales < 1 GeV, 
00 ! lattice simulations are the only tool to eventually give reliable answers, provided 

' that systematic errors are controlled. As we shall see, at present it is still a long 

way to achieve this goal. In fact, lattice investigations at finite density are hampered 
^ \ by the "sign problem", and only approximate methods are available that work at 

^ • small quark densities, ^ = ^b/'^'^T?^'^^ This adds further systematic errors to 

those known from zero density thermodynamics, like finite volume and discretisation 
effects. Accordingly, in this contribution we shall only consider the quark hadron 
transition at small densities. The widely accepted expectation is for a finite den- 
sity first order phase transition terminating in a critical endpoint, and an analytic 
crossover behaviour at /U = (cf. Fig. 3 (right)). 

In the absence of first principles calculations, where did this picture come from? 
It is based on combining lattice results for = in the larger parameter space 
{f^u,d-,'^siT} with model calculations at 7^ 0^^ and connecting various limiting 
cases by universality and continuity arguments.^) The schematic situation is depicted 
in Fig. 1. In the limits of zero and infinite quark masses (lower left and upper right 
corners), order parameters corresponding to the breaking of a global symmetry can 
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Fig. 1. Left: Schematic phase transition behaviour of Nf = 2 + 1 QCD for different choices of quark 
masses {mu,d,'rns) at fi = 0. Middle/Right: Critical surface swept by the chiral critical line as 
is turned on. Depending on the curvature, a QCD chiral critical point is present or absent. For 
heavy quarks the curvature has been determined*' and the first order region shrinks with fi. 



be defined, and one numerically finds first order phase transitions at small and large 
quark masses at some finite temperatures Tc(m). On the other hand, one observes 
an analytic crossover at intermediate quark masses, with second order boundary 
lines separating these regions. Both lines have been shown to belong to the Z{2) 
universality class of the 3d Ising model. ^^"^^ 

The "derivations" of the generally expected QCD phase diagram^^ state two 
crucial assumptions: a) the chiral transition for Nf = 2 is second order and thus 
in the 0(4) universality class, which implies the existence of a tricritical point at 
some strange quark mass m}p^; b) by switching on ^u, this point will continuously 
move to larger m<j until it ends up as a tricritical point for the Nf = 2 theory at 
some finite ^Mtric- Similarly, for small but non-zero rra^ ,^, the Z{2) chiral critical 
line would continously shift with /i until it passes through the physical point at ^u^;, 
corresponding to the endpoint of the QCD phase diagram. This is depicted in Fig. 1 
(middle), where the critical point is part of the chiral critical surface. Note, however, 
that there is no a priori reason for this. In principle it is also possible for the chiral 
critical surface to bend towards smaller quark masses, cf. Fig. 1 (right), in which case 
there would be no chiral critical point or phase transition at moderate densities. In 
the sequel the lattice evidence for these scenarios will be reviewed. 



§2. Nf = 2 at zero density 



Let us first consider assumption a) above. Since the cost of dynamical simu- 
lations explodes with shrinking quark masses, the order of the two-flavour chiral 
transition is even harder to determine than that of physical QCD. Despite many at- 
tempts it could not yet be conclusively settled. Wilson fermions appear to see 0(4) 
scaling,^) while staggered actions are inconsistent with 0(4) and 0(2) (for the discre- 
tised theory ).-'^'^) A recent finite size scaling analysis using staggered fermions with 
unprecedented lattice sizes was performed in.^^^ Again, these data appear inconsis- 
tent with 0(4)/0(2), and the authors conclude a first order transition to be more 
likely. A different conclusion was reached in,^^) in which xQCD was investigated 
numerically. This is a staggered action modified by an irrelevant term (i.e. one going 
to zero in the continuum limit) such as to allow simulations in the chiral limit. The 
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Fig. 2. Left: The chiral critical line in the bare mass plane at /i = 0. Nf = 3 corresponds to the 
solid line. Also shown is a fit to an assumed tricritical point, m*'"*'^ ~ 2.8r.^^-' Right: Finite 
size scaling of the chiral susceptibility at the physical point for Nt = 4, 6. The peaks saturate 
at a finite value, consistent with crossover behaviour. 



authors find their data compatible with those of an 0(2) spin model on moderate to 
small volumes. They thus suspect that finite volume effects of current Nj = 2 QCD 
simulations mask the correct scaling. 

Finally, from universality of chiral models it is known that the order of the chiral 
transition is related to the strength of the C/a(1) anomaly. ^'^^ In a model constructed 
to have the right symmetry with a tunable anomaly strength, it has recently been 
demonstrated non-perturbatively that both scenarios are possible, with a strong 
anomaly required for the chiral phase transition to be second order. 

Thus, we cannot yet take assumption a) for granted. Should the chiral transi- 
tion turn out to be first order, the likely modification of Fig. 1 (left) would be the 
disappearance of the tricritical point, with the chiral critical line intersecting the 
Nf = 2 axis at some finite m-u.d and being Z(2) all the way. Provided this line shifts 
to the right with as in assumption b), the (T, /i^)-phase diagram for physical QCD 
would still look as expected. But contrary to the scenario in Ref.,^) its critical point 
would be unrelated to any tricritical point. 

§3. The chiral critical line at = 

The boundary line between the chiral first order and crossover regions has re- 
cently been mapped out on A^^ = 4 lattices, Fig. 2 (left). A convenient observable 
is the Binder cumulant Ba{X) = {{X - {X)Y) / {{X - {X)f)'^, with X = ij^. At the 
second order transition, i?4 takes the value 1.604 dictated by the ScZ Ising universality 
class. In agreement with expectations, the critical line steepens in approaching the 
chiral limit. Assuming a tricritical point on the m^-axis according to Fig. 1 (left), 
the critical line is in fact consistent with tricritical scaling with rn-„ ^f^^ and allows 
to estimate m*"'^ ~ 2.8Tc. Note however, that this estimate is extremely cut-off 
sensitive and would change considerably on a finer lattice. 

The most important question concerns the location of the physical point, which 
is marked by the cross in Fig. 2 (left). As expected, it is on the crossover side of the 
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Left: Comparison of different methods to compute the critical couplings. Right: Phase 
fram for physical quark masses as predicted by two parameter reweighting on Nt — 4.^^-* 



critical line. In Ref. 15) ratios of pion to rho and kaon masses were evaluated on the 
points marked by arrows, to ensure that this statement indeed carries over from bare 
quark masses to the spectrum of physical particles. Quark masses being extremely 
susceptible to renormalisation effects, it is important to check this situation on finer 
lattices. This has been completed in Ref. 16) by a slightly different strategy. Here, 
the authors tune the quark masses to the physical theory and then perform a finite 
size scaling analysis of susceptibilities around the critical temperature, viz. lattice 
coupling. This is shown in Fig. 2 (right) for different lattice spacings. Clearly, the 
peaks saturate at a finite value which can be extrapolated to the continuum. This 
quite convincingly shows physical QCD to exhibit an analytic crossover rather than 
a true phase transition at zero density. The only remaining caveat to this conclusion 
would be if there was a fundamental problem with the so-called rooting trick when 
using staggered fermions, as frequently debated. ^''^ Similar calculations using Wilson 
fermions could close this gap soon. 

In an attempt to control cut-off effects, also the critical line has been checked on 
finer Nt = 6 (a ~ 0.2 fm) lattices. The results show an important shift of the critical 
line towards the origin: for the Nf = 3 theory, the pion mass, measured at T = 
with the critical quark mass, decreases from 1.6 Tc to 0.95 Tc.^^-* Similar results 
are reported for Nf = 2 + 1.^^-* Since cut-off effects on the physical point are much 
milder, this considerably increases the distance of the critical surface to the physical 
point, Fig. 1. Regardless of the sign of the curvature, this trend alone makes a QCD 
chiral critical point at small fi /T < 1 less likely. 



4. Calculations at finite density 



Straightforward Monte Carlo simulations at finite baryon density are impossi- 
ble. This is because the fermion determinant becomes complex for non-vanishing fj,, 
prohibiting its use as a probability weight in Monte Carlo algorithms. This fact is 
also known as the "sign-problem". 

There is a number of methods that circumvent the sign problem, rather than 
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Table I. Coefficient t2 in the Taylor expansion of the transition line, Eq. (4-1) from Nt = 4. 



solving it: i) Multi-parameter reweighting, ii) Taylor expansion in {/jl/T)'^ around n = 
(even powers because of CP invariance) and iii) simulations at imaginary chemical 
potential, either followed by analytic continuation or Fourier transformed to the 
canonical ensemble. All of these introduce some degree of approximation. However, 
the systematic errors are rather different, thus allowing for powerful crosschecks. 
Reviews specialized on the technical aspects can be found in Refs. 2), 3). 

The first task is to identify the phase boundary, i.e. the critical coupling and thus 
Tc{n). This has been done for a variety of flavours and quark masses using different 
methods. For a quantitative comparison one needs data at one fixed parameter set. 
Such a comparison is shown for the critical coupling in Fig. 3 (left), for Nf = 4 
staggered quarks with the same action and quark mass m/T = 0.2. (For that quark 
mass the transition is first order along the entire curve). One observes quantitative 
agreement up to fi/T ~ 1.3, after which the different results start to scatter. Thus, 
all methods appear to be reliable for n/T <1, or fis 500 MeV. The case of physical 
quark masses, after conversion to continuum units, is shown in Fig. 3 (right). One 
observes that Tc is decreasing only very slowly with This is consistent with a 
description by a series in (/x/ttT)^ with coefficients of order one, 

Z|).,_,(.„„,,(A)Vo((A)') . ,4.1) 

The leading coefficients for various cases have been collected from the literature^^ 
and are reproduced in Table I. The curvature grows with Nf, which is consistent 
with ~ Nf/Nc behaviour found in large Nc expansions. Subleading coefficients 
are emerging at present but not statistically significant yet. Note that continuum 
conversions relying on the two-loop beta function are certainly not reliable for these 
coarse lattices, while fits to non-perturbative beta functions tend to increase the 
curvature. 



§5. The chiral critical surface 



All methods mentioned here also give signals for criticality, but the comparison 
is non-trivial because of different parameter sets. A simulation using reweighting 
methods on A/^ = 4 lattices puts the critical point at 11% ~ 360 MeV,23) Fig. 3 
(right), supporting the standard expected scenario. Quark masses were tuned to 
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Fig. 4. Left: Chiral critical line at zero and non-zero imaginary /i. Right: One sigma error band 
for the Nf = 3 critical mass am'^{fi) resulting from a leading order fit. Both for Nt = 4.^^-' 



give the ratios niT^/mp ss 0.19, m7r/mi<: ^ 0.27, which are close to their physical 
values. In principle the determination of a critical point is also possible via the 
Taylor expansion, where a true phase transition will be signalled by a finite radius 
of convergence for the pressure series about = as the volume is increased. A 
critical endpoint for the Nj = 2 theory, based on this approach, was reported in^^^ 
for bare quark mass m/Tc = 0.1. Taking the measured first four coefficients for the 
asymptotic behaviour of the series, the estimate for the location of the critical point 
is ii^/Te = 1.1 ± 0.2 at Te/Tc{ij- = 0) = 0.95. 

Rather than fixing a theory with a particular set of quark masses and then 
switching on a chemical potential, let us now try to learn about the phase structure 
in the extended parameter space {m^^^, m^, T, /i}, i.e. map out the chiral critical 
surface. This has been done for Nt = 4 lattices using simulations at imaginary 
chemical potential. Fig. 4 shows a comparison of the chiral critical line at zero 
density and a few points at imaginary chemical potential = i2AT. The finite 
density effect is very small, consistent with what is found for the change of the 
critical temperature, Eq. (4-1). Thus, the critical surface in Fig. 1 appears to emerge 
very steeply from the quark mass plane, making the critical point of physical QCD 
extremely quark mass sensitive. 

The critical surface for imaginary /x. Fig. 4, is moving to larger quark masses. 
What does this imply for real chemical potential? To answer this question we focus 
on Nf = 3, collecting data for several values of /ij. Since ^Uj/T < 1, the critical quark 
mass may be Taylor expanded anf{^) = om^ + c'^(a/x)^ + . . ., and the coefficients 
can be fitted to the data at imaginary /x. Fig. 4 (right) shows a one sigma error 
band for the critical bare quark mass from a leading order fit. As observed before, 
the //-dependence is very weak and even consistent with zero for these errors. Since 
T{^) = l/{a{fx)Nt), a{fj,) is an increasing function on a given lattice, so that the 
critical mass in fixed physical units shrinks with //, as in Fig. 1 (right). 

One may worry about systematic errors when fitting a leading order poly- 
nomial to data containing the full functional dependence. For example, subse- 
quent terms may be cancelling in the imaginary, but not in the real direction. 
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Fig. 5. Leading /^^-derivatives of the critical quark mass on A't = 4 on L = 8 (left) and L — 12 
(right). The error bands give the corresponding values of LO and NLO fits to Taylor series. 



In order to check this, we have also calculated the leading derivative directly via 
di = —d(ai_i)2B4^/damB4. We do this in a novel, efficient way by evaluating finite 
differences 

Because the required shift in the couplings is very small, it is adequate and safe to 
use the original Monte Carlo ensemble for amg, ^ = and reweight the results by the 
standard Ferrenberg-Swendsen method. Moreover, by reweighting to imaginary fi the 
reweighting factors remain real positive and close to 1. The results of this procedure 
based on 5 million trajectories on 8^, 12^ x 4, are shown in Fig. 5 for two volumes. 
Subleading terms show up as a slope in the linear extrapolation. Indeed, such a slope 
is visible in Fig. 5 (right). Both coefficients are consistent with the results from the 
finite /ij calculations, provided the next-to-leading order is taken into account. We 
have also continued to collect statistics for the imaginary /x simulations, so we now 
have two significant terms. Putting everything together, our current best estimate 
for Nf = 3 on iVj = 4 lattices is 

!^£M = (-JLf _„m (JlV 



The emerging ^^-term is negative as well, further shrinking the first order region with 
increasing /i. On ./Vf = 4 lattices, there is thus little doubt that the scenario Fig. 1 
(right) is realised for ^b<500 MeV. Note that this is analogous to the situation 
with heavy quarks, where the first order transition is also weakening with /i.^-* 

§6. Discussion 

Our result for the chiral critical surface appears to be in contradiction with the 
phase diagram obtained from reweighting methods. Fig. 3. As explained above, we 
have checked with independent methods that the signs we find for our coefficients 
are not artefacts of the fitting procedure. On the other hand, there are concerns 
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that the critical point determined via reweighting fahs into a parameter range where 
reweighting becomes problematic.^^) However, even disregarding those, the bare 
lattice results of Rcfs. 15), 23) are not necessarily inconsistent. Reweighting in fi 
is performed for fixed quark masses in lattice units, aniq. Since T(^) = l/{aNt), 
a{iJi) is an increasing function on a given lattice, hence the critical point observed 
by reweighting corresponds to quark masses smaller than physical. Indeed, Fig. 2 
shows that on iVt = 4 the physical point is very close to the critical surface. This 
is a discretisation effect, and simulations on finer lattices are required in order to 
settle the issue. As discussed above, we know already that the distance between the 
critical surface and the physical point grows significantly, making a chiral critical 
point at small hb less likely, irrespective of the curvature of mc{lJ-)- Calculations of 
the curvature on Nt = 6 are currently in progress. 
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